Nuclear transparency from quasielastic 12 C(e,e'p) 
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We studied the reaction 12 C(e,e'p) in quasielastic kinematics at momentum transfers between 0.6 
and 1.8 (GeV/c) 2 covering the single-particle region. From this the nuclear transparency factors are 
extracted using two methods. The results are compared to theoretical predictions obtained using a 
generalization of Glauber theory described in this paper. Furthermore, the momentum distribution 
in the region of the ls-state up to momenta of 300 MeV/c is obtained from the data and compared 
to the Correlated Basis Function theory and the Independent-Particle Shell model. 

PACS numbers: 21.10.-k, 25.30.Dh 



I. INTRODUCTION 

Electromagnetic probes are a valuable tool to examine 
the structure of nuclei. The distortion of the electron is 
small and the entire nucleus is probed. In (e,e'p) reac- 
tions one can assume in good approximation that only 
one proton is involved in the primary reaction. However, 
on the way out of the nucleus the nucleon is subject to 
reactions which can remove it from the reaction chan- 
nel under consideration. Such inelastic processes lead to 
absorption and deflection of the outgoing proton and con- 
sequently to a reduction of the detected yield. At high 
outgoing proton energy, the main process of concern is 
the reduction of the proton flux. The reduction factor is 
called nuclear transparency Ta- An understanding of the 
propagation of nuclcons in nuclear matter is important 
for the interpretation of many experiments. 

There are two theoretical approaches to deal with ab- 
sorption. For low proton kinetic energies (T p < 200 MeV) 
an optical potential is used. The parameters of the poten- 
tial are chosen such as to reproduce the phase shifts and 
cross sections of nucleon-nucleus scattering. Rs imagi- 
nary part accounts for inelastic processes, i.e. absorp- 
tion. For nucleon momenta larger than 1 GeV/c where 
the inelastic part of the free nucleon-nucleon (NN) cross 
section dominates, Glauber calculations are often used. 
They derive the Final State Interaction (FSI) of the prop- 
agating nucleon in nuclear medium directly from the el- 
ementary nucleon-nucleon cross section cat at. In the nu- 
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clear medium the N-N cross section has a similar energy 
dependence as the free one for T p larger than 1 GeV 
but is reduced by about 20%. This value is roughly in 
agreement with the one obtained using a simplified geo- 
metrical model assumingclassical attenuation and fitting 
it to experimental data []j . At low T p stronger modifica- 
tions of the N-N cross section occur due to Pauli Blocking 
and dispersion effects as it was shown in Q. Therefore 
Glauber calculation are generally thought to have a lower 
limit of validity in the nucleon kinetic energy. 

In this paper it will be shown that agreement with the 
data down to T p = 0.3 GeV can be achieved using an ap- 
proach based on the correlated Glauber approximation 
to take into account these effects. A similar approach 
was already able to describe the Ta measurement at 
180 MeV for 12 C, 27 Al, 58 Ni and 181 Ta within the exper- 
imental error bars. Details of the theoretical calculation 
are given in section [H] The theoretical results arc com- 
pared to data taken as part of the experiment E97-006 
and previously measured values for the nuclear trans- 
parency described in section 1 1 1 1 1 For the extraction of 
Ta an improved method was applied in this paper and 
the transparencies obtained via the standard analysis are 
given for comparison. The improvement is not only ap- 
parent in the better agreement of Ta with the theory but 
also when comparing the momentum distribution in the 
region of the ls^-statc in 12 C with the Independent- 
Particle Shell Model (IPSM) and the Correlated Basis 
Function theory (CBF) For this deeper lying state 
the deviation from the shape of the IPSM momentum dis- 
tribution is already apparent at momenta > 150 MeV/c. 
However, the CBF theory gives a good description up to 
250 MeV/c. 
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II. THEORY 



Neglecting many-body contributions to the target elec- 
tromagnetic current, the nuclear matrix element entering 
the definition of the (e, e'p) transition amplitude can be 
written (see, e.g., Rcf. |f|) 



M„(p,q) = (*< 



E4- 



q ak|*o) 



(1) 



where p and q denote the momentum of the detected 
proton and the momentum transfer, respectively. In the 
above equation, and are nucleon creation and an- 
nihilation operators, \^f ) is the target ground state, sat- 
isfying the many-body Schrodinger equation 



H A |*0> = £*>|*0> , 



(2) 



and is the final scattering state. 

The effect of FSI between the knocked out nucleon 
(labelled with index 1) and the spectator particles can be 
best analyzed rewriting the A-body nuclear hamiltonian 
Ha in the form 



H A — H + H FS i — (-Ha-i + Ti) 



H 



FSI 



(3) 



where i?A-i and T\ are the hamiltonian of the recoiling 
(A-l)-nucleon system and the kinetic energy operator 
associated with the struck nucleon, respectively, while 
i^FSi describes the interactions between the struck nu- 
cleon and the spectators. Obviously, PWIA amounts to 
setting TJpsi = 0. 

The decomposition of Eq.|j3J) can be used to write 
\^kp) in the form [(| 



|^p ) ) = njr ) l*np> , 



(4) 



where the asymptotic state |$ n p)j describing the system 
in absence of FSI, is an eigenstate of Hq of the product 
form 



l$np) = |P) ® \<Pn) , 



(5) 



p) and \(f n ) being eigenstates of T\ and Ha-i, respec- 
tively. 

In coordinate space l^np) can be written (in order to 
simplify the notation spin indices will be omitted) 



(6) 



(R\<S>n P ) = $np(R) = \ TJ e'P >„(]i) 



where R = {ri, r 2 . . . , r^} and R = {r 2 , r^} spec- 
ify the configurations of the full A-particle system and the 
(A-l)-particle spectator system, respectively, whereas V 
denotes the normalization volume. 

Setting Hp ^ = 1, which amounts to disregarding the 
effects of FSI, and substituting the resulting final state 
into Eq.QJ, one obtains the PWIA transition amplitude, 
that depends upon the initial momentum k = p — q only. 



The scattering operator fl p ' in EqQJdcscribes the dis- 
tortion of the asymptotic wave function resulting from 
rescattering of the knocked-out nucleon. It can be for- 
mally written as [|J 

fi<T> = lim e tHAt e- %Hat 

p t^oo 



= lim T c 

t — >00 



-i f*dt' Hpsi(f') 



(7) 



where T is the time ordering operator and 

H FSI (t) = e iH ° t H FSl e- iH ° t . (8) 

The calculation of Q p ' from Eq.Q using a real- 
istic nuclear hamiltonian involves prohibitive difficul- 
ties. However, when the kinetic energy carried by the 
knocked-out proton is large, the structure of fl p can 
be greatly simplified using a generalization of the ap- 
proximation scheme originally developed by Glauber to 
describe proton-nucleus scattering Q. 

The basic assumptions underlying this scheme, gen- 
erally referred to as correlated Glauber approximation 
(CGA) (see, e.g., Ref. @ and references therein), are 
that i) the fast struck nucleon moves along a straight 
trajectory, being undeflected by rescattering processes 
(eikonal approximation) and ii) the spectator system can 
be approximated by a collection of fixed scattering cen- 
ters (frozen approximation). 

Within CGA, Q p ~^ can be written in coordinate space 



as 



(R\n^\R) = n p -\R) = p z 



i-£r P (i,j) 

J=2 



r p (i,j)T p (i,k) 

k>j=2 



(9) 



where the positive z-axis is chosen along the cikonal tra- 
jectory and the z-ordering operator P z prevents the oc- 
currence of backward scattering of the fast struck proton. 
The profile function r p , appearing in Eq.@ is defined 

as 



r p (i,j) = e(z j -z 1 ) 7p (\b 1 -b j \) 



(10) 



where the step function preserves causality and |bi — hj\ 
is the projection of the interparticle distance on the im- 
pact parameter plane (i.e. the xy plane). The function 
7 p (fe), containing all the information on the dynamics of 
the scattering process, is simply related to the NN scat- 
tering amplitude at incident momentum p and perpen- 
dicular momentum transfer k t , f p (k t ), through 



d 2 k 



(2nf 



- e ikt b 



f P (h 



(11) 



At large p, the scattering amplitude f p (kt) ex- 
tracted from the measured NN cross section is usually 
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parametrized in the form 



f P (h) = i o- pN (l - ia pN )t 



(12) 



where <j p n and a P N denote the total cross section and 
the ratio between the real and the imaginary part of the 
amplitude, respectively, while the slope parameter B is 
related to the range of the interaction. In case of zero- 
range, corresponding to l/B = 0, the impact param- 
eter dependence of Jp(b) reduces to a two-dimensional 
(^-function. 

The calculation of the scattering operator of Eq. (9) 
requires the knowledge of the NN scattering amplitude 
in the nuclear medium. Pandharipande and Pieper Q 
have shown that the total NN cross sections in nuclear 
matter and in vacuum can be easily related to one an- 
other. Their approach, in which the velocity dependence 
of the nuclear mean field and Pauli blocking of the final 
states available to the spectator nucleons are both taken 
into account, allows one to obtain the effective total cross 
section, a P N < <j p n , from the measured differential cross 
section da p ^/dD,. The results of Ref. @ suggest that 
medium effect reduce the total cross section by as much 
as ~ 50 % and ~ 18 % at incident energies of .2 and 1. 
GeV, respectively. 

Theoretical studies of the (e, e'p) cross section carried 
out using CGA lCj show that at large proton momentum 
the dominant FSI effect is a quenching that can be ac- 
curately described introducing a transparency factor Ta, 
written in terms of the scattering operator of Eq.© as 



5] 



d\p p (v) \Q ( p >(r 



(13) 



In the above equation p p (r) is the proton density of the 
target, normalized to Z, and 



Pp( r i)l^p (ri) 



dR |*o(fl)| a W ■ (14) 



The full calculation of the right hand side of Ea. (|14|) . 
involves 3(A-l)-dimcnsional integrations. It has been 
carried out for 4 He and 16 using Monte Carlo techniques 
and realistic many-body wave functions 0, llOj . 

Approximations to Ta can be obtained expanding the 
integrand of Eq. (|T3")l according to 



1 



ft,(ri)|r£'(ri)| a 
1 



J rf 3 r 2 r(l,2)^(r l5 r 2 ) 



Pp(ri) 

- j d 3 r 2 d 3 r 3 r(l,2)r(l,3)^(r 1 ,r 2 ,r 3 ) 



+ / d 3 r 2 d 3 r 3 d 3 r4 



(15) 



where the distribution functions p p n ^ ^(ri, r 2 . . . r„) 
yield the joint probability of finding the struck proton 
at ri and the n — 1 spectator nucleons at r 2 . . . r„. 



The quantity in square brackets in Ea. (|15f) describes 
FSI effects, that lead to a departure from the PWIA re- 
sult Ta = 1. The contribution of processes in which the 
struck proton undergoes n — 1 rescatterings involves the 
7i-nucleon distribution. For example, the single rescatter- 
ing term involves the two-nucleon distribution that can 
be written in the form j 1 1 1 | 



Pp 2 i(ri,r 2 ) = p p (r 1 )p N (r 2 )g(r 1 ,r 2 ) 



(16) 



The function <7(ri,r 2 ) describes the effects of dynamical 
correlations. Due to the strongly repulsive nature of NN 
interactions at short range, g(r) <C 1 at r = ri — r 2 | <^ 1 
fm, while at large r g(r) — > 1 and the distribution func- 
tion reduces to the prediction of the independent particle 
model. 

Equation (|15|) shows that inclusion of NN correlations 
produces a nontrivial pattern of effects. For example, the 
presence of g(ri, r 2 ) in the two-body (single rescattering) 
term enhances the transparency, since the short range 
repulsion between the struck particle and the spectator 
reduces the rescattering probability. On the other hand, 
the repulsion between two spectators in the three-body 
(double rescattering) term leads to the opposite effect. 

As short range correlation are known to be largely un- 
affected by shell and surface effects, <?(ri,r 2 ) can be ob- 
tained using the Local Density Approximation (LDA), 
i.e. setting 



ff(ri,r 2 ) « 9nm 



\vi - r 2 | , p N 



ri + r 2 



(17) 



where gNAi(r 7 p) is the radial distribution function calcu- 
lated in uniform nuclear matter at constant density p. 

A similar procedure can be employed to construct 
distribution functions involving three or more nucleons, 
needed to evaluate the contributions of processes involv- 
ing more than one rescattering. The LDA scheme, com- 
bining nuclear matter results and the measured one-body 
density, provides a consistent framework to evaluate Ta 
from Eqs. (|13|l and l|15(l provided the NN scattering am- 
plitude is known. 

The theoretical curves shown in Fig. 4 have been ob- 
tained using LDA and an effective NN scattering ampli- 
tude written as in Eq. I|12|) . Medium effects have been 
taken into account replacing the total cross sections in 
vacuum with those of Ref. 0, while for the slope pa- 
rameters the free space values resulting from the fit of 
Ref. ^3] j based on the data of Ref. 01 j have been used. 
The third parameter appearing in Eq. I|12|). namely a P N, 
does not affect the transparency factor. 

The calculations have been carried out including two- 
, three- and four-body terms in Ea. (|15fl . corresponding 
to single, double and triple rescattering. The inclusion 
of three- and four-body contributions to Ta produces a 
change of ~ 10% and ~ 3%, respectively, in 12 C and 
~ 15% and ~ 5% in 197 Au. 
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TABLE I: For the five settings in quasielastic kinematics the 
averaged momentum transfer Q 2 , the beam energy E e , the 
central momentum of the HMS p e i (electron) and of the SOS 
p p (proton) as well as the central angles 9 e and 6 P for HMS 
and SOS are given. 



Q 2 


E e 




0, 


P P 


Op 


(GeV/c) 2 


(GeV) 


(GeV/c) 


degree 


(GeV/c) 


degree 


0.59 


3.298 


2.95 


14.4 


0.85 


60.3 


0.80 


3.298 


2.75 


17.0 


1.00 


56.2 


1.13 


3.123 


2.50 


22.2 


1.25 


49.7 


1.53 


3.298 


2.40 


25.4 


1.50 


44.6 


1.85 


3.298 


2.28 


29.0 


1.70 


40.7 



III. EXPERIMENT AND DATA ANALYSIS 

The experiment was performed at the Thomas Jeffer- 
son National Accelerator Facility (TJNAF) in hall C. 
Electrons were detected in the High Momentum Spec- 
trometer HMS, protons in the Short Orbit Spectrometer 
SOS. Data were taken in five kinematics shown in Tab. [fl 
on a 12 C target of 2.5% radiation length thickness. The 
target thickness was determined by measuring the weight 
and the size of the target. Its uncertainty is estimated to 
be 0.3%. The data cover a range of Q 2 corresponding to 
T p = 0.3 to 1.0 GeV. In the same kinematics data with a 
liquid hydrogen target were taken. This target consists of 
an upright standing aluminum cylinder with a diameter 
of 4 cm and 0.13 mm thick windows. The background 
contribution from the windows was subtracted using two 
aluminum targets spaced by the dimension of the cryo- 
genic target. These data serve as a check of the analysis 
and for the determination of the kinematical offsets to 
the values given in Tab. [I] The offset to the beam energy 
was determined to -0.1% from an independent analysis 
of ^3 using a large set of H(e,e'p) data. From the re- 
measurement of the magnetic dipole field in the magnets 
located in the arc of the beam line to hall C which are 
used for the determination of the beam energy, a correc- 
tion of -0.2% was found Both offsets were used in 
the following analysis and contribute to the systematic 
error. 

The H(c,c'p) data were also used to determine the frac- 
tion of protons absorbed in the target and on their way 
to the detector system. For this the rates of the reac- 
tion H(e,e'p) and H(e,e') in the HMS were compared in a 
kinematical region covered by both reactions. The result 
of (0.953 ± 0.011) agrees with the one calculated from 
the mean free path and the material thicknesses. 

The analysis of the H(e,e'p) data revealed that for mo- 
menta larger than 1 GeV/c the protons punch through 
the collimator which consists of 2.5 inch tungsten. This 
leads to additional energy loss of 50 to 300 MeV for less 
than 3% of the events. The effect was included in the 



simulation (see below). 

The H(e,e'p) data are in good agreement with previous 
results. The momentum acceptance relative to the cen- 
tral momentum Ap/p was -8.8% to 10.4% for the HMS 
and ± 15% for the SOS. It was determined using overlap- 
ping spectra from inclusive electron scattering on 12 C in 
the inelastic regime where the cross section is a smooth 
function of the momentum. Both spectrometers were 
fixed at 20° and the central momentum was varied in 
steps of 0.1 GeV/c. This guarantees an overlap between 
the settings. A correction of the order of ± 2% was ap- 
plied to the yield as a function of Ap/p. 

The spectra from the H(e,e'p) data were compared to 
the Monte Carlo simulation SIMC of the hall C collabo- 
ration and good agreement were found. The simulation 
describes the electron beam entering the target (includ- 
ing the rastering of the beam over the target used to avoid 
local heating) and follows the particles from the reaction 
vertex to the detector system. It employs transfer ma- 
trices for the spectrometers and takes energy loss and 
multiple scattering in material into account. Radiative 
corrections are calculated according to ^(|. The effect 
of Coulomb distortion on the electron wave function in 
the vicinity of heavy nuclei is small at the momenta con- 
sidered here (~ 1%). It was taken into account by using 
the effective momentum approximation. The angular and 
momentum resolutions were adjusted to reproduce best 
the measured spectra. For this the spectra of the miss- 
ing energy E m and missing momentum p m are suitable 
in particular. A resolution (FWHM) in E m of w 7 MeV 
and in p m of « 9 MeV/c is achieved. From the measured 
kinematic variables, the missing energy is reconstructed 
according to 



EL - T n - T, 



R- 



(18) 



Here, E e is the beam energy and E e i and T p are the 
(kinetic) energies of the outgoing electron and proton, 
respectively Tr is the kinetic energy of the (undetected) 
recoiling (A-l)-system, which is reconstructed from p m 
using the spectator model. The missing energy can be 
identified with the removal energy E of the nucleoli. In 
the single-particle region it is the energy needed to re- 
move the nucleon from a particular state within the nu- 
cleus. In Plane Wave Impuls Approximation (PWIA) the 
initial momentum k of the nucleon bound in the nucleus 
is equal to the opposite of the missing momentum p m = 
q — p where q is the momentum transfer from the elec- 
tron to the nucleon. In PWIA the 6-fold differential cross 
section 



d 6 <r 



dE dE' dn e dE p dfl p 



K a ep S(E, k) Ta{Q 2 ) (19) 



factorizes into an elementary e-p cross section <j ep and 
the spectral function S(E, k) containing the information 
about the nuclear structure. The spectral function is the 
probability to find in the nucleus a nucleon with energy 
E and momentum k. 
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FIG. 1: Comparison of the experimental yield obtained as 
a function of E m in the single-particle region to the Monte 
Carlo simulation (solid) for the kinematics at Q 2 — 1.13 
(GeV/c) 2 (s. Tab. QJ. Both contributions contain radiative 
processes. 



The choice of cross section <r ep for the scattering of an 
electron on a bound nucleon is not unique. Since the 
off-shell form factors are not known the same hadronic 
currents as in the case of free nucleons are used but do not 
lead to the same result. The deviation between different 
choices increases with E and k. For the result of the 
present work, which is restricted to the single-particle 
region, the ambiguity is small. To be consistent with 
previous analyses the off-shell prescription er^ 1 of ref. 
[13 was used. The deviation to af p 2 [lj| and to af p 
is less than 1.5% for the five settings in the E m -p m region 
used for this analysis. 

As parameterization for the electromagnetic form fac- 
tors of the proton we use the dispersion-theoretical anal- 
ysis based on vector meson dominance model of |l8j . 
The change on the cross section at the highest Q 2 of 
the present work is less than 0.5 % when the result from 
the double-polarization method for the G ep /G mp ratio 
0,123] is employed. 

In the framework of PWIA the nuclear transparency 
enters as correction factor on the right hand side of 
Eq. (|19fl . Experimentally it is obtained by comparing 
the measured yield N exp (E m , p rn ) to the simulation 
N sm (E m , p m ) integrated over the same phase space V 



T A (Q 2 



J v dp m dE m N ex P(E m ,p m ) 
Ldp m dE m N sim (E m ,p m ) i 



(20) 



To restrict the analysis to the single-particle region cuts 
in E m < 0.08 GeV and p m < 0.3 GeV/c were applied. 
The same cuts were also used in most of the previous 
experiments. 

In the simulation a spectral function has to be chosen 



which correctly describes the experimental distribution 
well in the region of interest. In Fig. ^the i? m -spectrum 
of the experiment (dots) is compared to the simulation 
using an IPSM spectral function as input (solid curve). 
This spectral function factorizes into an E and k distri- 
bution for each orbit. It was also used for the analyses 
of the experiment NE-18 |2~3] at SLAC and experiments 
at TJNAF 0, El]. The ^-distribution is described by 
a Lorentzian and the momentum distribution is derived 
from a Wood-Saxon potential whose parameters are ad- 
justed to data measured at Saclay |24|. For this figure 
Ta = 0.6 in the simulation was used. The difference in 
the yield reflects the depletion of the single-particle re- 
gion due to short-range correlations (SRC). To account 
for this the simulated yield for carbon is divided by a fac- 
tor of e SRC 1.11 ± 0.03. This factor was used in all the 
previous analyses and was estimated from an early exam- 
ination of N-N correlations in 12 C and 16 [H|26|]. The 
inverse of e SRC corresponds to the occupation number, 
which modern many-body theories predict to be lower, 
« 80% HH3. This number is larger than the spectro- 
scopic factor giving the occupation of a state; it contains 
the background term induced by correlations which can- 
not be attributed to a specific orbit [2l| ■ The above 
considerations emphasize that the number of nucleons 
missing due to absorption in the nucleus cannot be dis- 
tinguished from the depletion of the single-particle orbits 
due to SRC. 
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FIG. 2: _B m -distribution for a slice in p m of 250 ± 50 MeV/c 
in the quasi-elastic kinematics at Q 2 = 1.53 (GeV/c) 2 . The 
experimental yield is shown as data points, the other curves 
are simulated taking radiative processes into account. Using 
the IPSM model results in the lower curve (dashed) . For the 
distribution (solid) the spectral function of Q is used as input 
for the Monte Carlo. 



To improve upon this approach a spectral function con- 
taining SRC from the beginning is employed in Ea. (|20ll . 
This spectral function is composed of a part due to SRC 
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0.05 0.1 0.15 0.2 0.25 0.3 
P m (GeV/c) 

FIG. 3: Momentum distribution in the region of 0.03 GeV < 
Em < 0.08 GeV obtained from the data taken in the kinemat- 
ics of Tab. □ (squares), the CBF theory (solid) and the IPSM 
(dashed). Three data poin ts (circles) are from data focusing 
on the high p m region |2l], [33 • 



which accounts for 22% of the total strength as calculated 
in LDA, and the IPSM spectral function mentioned above 
but reduced by a factor (1 - 0.22) to ensure normaliza- 
tion. This approach circumvents the application of the 
extra factor e which is a poor approximation because 
the amount of strength due to SRC is increasing with p m . 
Its limitation is demonstrated in Fig. [3 For a p m bin 
of (250 ± 50) MeV/c the experimental -E m -distribution 
is compared to the simulation using the IPSM spectral 
function (including the correction e SRC ) and the CBF 
spectral function. It is obvious that the IPSM fails al- 
ready at moderate p m and E m whereas the CBF theory is 
in good agreement with the data. Note that the increase 
of the yield at larger E m is not due to S(E, k), but is due 
to bremsstrahlung shifting events from low to high E m . 

IV. RESULTS 

When extracting Ta using Eq. good agreement 
between Monte Carlo simulation and data is required. 
Spectra of kinematical quantities like momenta and an- 
gles depend mainly on the capability to map the optics of 
the spectrometers. These were verified using the H(e,e'p) 
data. The comparison with the simulated E m and p m dis- 
tribution then reflects the quality of the spectral function 
used as input for the simulation. The agreement in the 
-Em distribution can easily be judged from the spectrum 
given in Fig. ^ In contrast to previous works, where the 
same spectral function was used, the width of the I7J3/2 
state had to be reduced from 5 MeV to 1 MeV. This 
value is closer to the expectation of a vanishing width for 
the IP3/2 state. The momentum distribution of the ls 1 / 2 
state was obtained from the data. For this the p m distri- 
bution was split into 10 MeV/c bins up to 300 McV/c. 
The data then were compared bin for bin to the simu- 
lation using the CBF spectral function. Simultaneously 



radiative corrections were taken into account by the de- 
radiation technique. The ratio was taken as a correction 
factor to the momentum distribution. The total normal- 
ization was adjusted to the simulation (using Ta = 0.6) 
which compensates for small differences in Ta in the five 
kinematics. The result obtained from the data (squares) 
is shown in Fig. [21 compared to the momentum distribu- 
tion from the CBF theory (solid blue line) and from the 
IPSM (dashed black line) integrated over the same re- 
gion in E m (0.03 - 0.08 GeV). The statistical error bars 
are smaller than the symbols. The systematic error is 
given below. Fig. [21 shows clearly that the IPSM fails 
already above p m = 0.15 GeV/c. This is in contrast to 
the measurements done at NIKHEF on various nuclei (s. 
refs. IlllEI) where good agreement was found with the 
momentum distribution derived from a Wood-Saxon po- 
tential up to momenta of 250 MeV. In these experiments 
only valence states at the Fermi edge were examined. 
For deeper lying states like the Is 1/2 state the influence 
of SRC leads to modifications of the shape already at 
moderate p m . 

Since no correction factor e is applied to the IPSM 
momentum distribution in figure [21 it exceeds the exper- 
imental momentum distribution by ss 20% at small p m . 
The deviation at large p m is not a too serious problem 
because most of the yield used for the extraction of Ta 
comes from smaller p m . The CBF theory gives a res- 
onable description to the data up to p m f=a 0.25 GeV/c 
which is close to the upper limit used in the Ta analysis. 
At higher p m the CBF theory underestimates the data in 
this E m region. 

The data points shown as full circles were also taken 
in the E97-006 experiment but with the aim to examine 
SRC at high E m and p m [2lll32T |. For this a spectral func- 
tion was extracted from the data. The momentum distri- 
bution was obtained by integrating the spectral function 
over the given E m region. These data are not subject 
of this article but the agreement in the overlap region 
of both data sets provides an important confirmation of 
consistency using different analysis methods. 

TABLE II: Nuclear transparency for carbon obtained using 
the IPSM and the CBF spectral function in the analysis. The 
error given is the quadratic sum of statistical and systematic 
uncertainty neglecting the model-dependent error. 



T p (GeV) T A (IPSM) T A (CBF) 

0.328 0.598 ± 0.023 0.641 ± 0.025 

0.433 0.581 ± 0.023 0.628 ± 0.025 

0.625 0.566 ± 0.022 0.605 ± 0.024 

0.830 0.561 ± 0.022 0.593 ± 0.023 

1.00 0.553 ± 0.022 0.591 ± 0.023 



The results for the nuclear transparency Ta obtained 
using the IPSM and the CBF spectral functions in 
Ea. (|20"|) are shown in table [B] When using the CBF 
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FIG. 4: Nuclear transparency Ta for C, Fe and Au as a func- 
tion of the proton kinetic energy T v compared to the corre- 
lated Glauber calculations (solid lines) . The data indicated by 
circles are from the NE18-experiment at SLAC [2^ |. squares 
and diamonds are Jlab data of l23ll and Q and from Bates 
Q (triangle down). The result indicated by stars is obtained 
with the correlated spectral function of 0. 



theory the results are about 5% higher than the results 
using IPSM. This can be traced back to the larger in- 
fluence of SRC predicted by the CBF theory, the corre- 
lated strength of which is in agreement with the strength 
found at high E and k by experiment |32l |. From the 
e SRC q UO ted above, which corresponds to a depletion in 
the single-particle state of 0.1, and the contribution of 
SRC in the CBF theory (0.22) one could naively have 
expected a larger difference between the two analyses. 
The deviation is not so large due to strength from SRC 
which also contributes to the single-particle region. 

It should be noted that the results for the momentum 
distribution and the nuclear transparency factor might 
depend slightly on the kinematics chosen. In general the 
yield obtained in perpendicular kinematics exceeds the 
one measured in parallel kinematics due to additional re- 
action mechanisms present in perpendicular kinematics 
[33j . In the single-particle region (below the Fermi mo- 
mentum) a difference of 6 % is expected [3 E3 ■ 

In the region of Q 2 = 0.15 - 0.6 (GeV/c) 2 refs. [H 
and found an excess of the transverse response 40 % 
in 12 C for E m > 0.025 GeV compared to a free proton. 
At higher Q 2 the transverse response quickly decreases. 
For the data point at Q 2 = 0.59 GeV/c) 2 this would 
lead to an increase of the transparency factor of w 9 %. 
To search for this effect in the data we extracted the 
transparency factor from the lp-state alone (E m < 0.025 
GeV). We find a decrease of Ta by 4 %, i.e. a two times 
smaller effect. At Q 2 = 0.8 (GeV/c) 2 the effect is 1 % 
and decreases futher for higher Q 2 . 

The error given in table [H] is the quadratic sum of 
the statistical and systematic uncertainties where the 
latter dominates with 3.9% in total. It consists of the 
already mentioned uncertainties in the target thickness 
(0.3%), the charge measurement (1%), proton transmis- 



sion (1.1%) as well as the uncertainties in several correc- 
tion factors that have to be applied to the experimental 
yield like dead time and detector efficiency (2%). The 
stability of the result against different cuts in the data 
was conservatively estimated to 2%. For the simulation, 
in particular the phase space, an error of 2% was taken 
into account. The influence of the choice of the kincmat- 
ical offsets was tested with the above mentioned second 
set of offsets corresponding to a difference in energy of 
0.1%. It leads to differences in the result of less than 
0.5%. For the fluctuations from run to run an uncer- 
tainty of 0.5% was estimated. The systematic error does 
not contain the model-dependent uncertainty. For the 
choice of the spectral function an error of 2% is taken 
into account which reflects mainly the agreement between 
data and simulation discussed above. The error in the 
correction factor e SRC of 3% has to be applied only in 
the case of the analysis using the IPSM spectral function. 
To estimate the influence of the off-shell cross section the 
analysis was repeated using the cc2- version of and the 
cc- version of [2l|,[32]. The latter leads to almost the same 
result as the ecl-version of 0. An error of 2% was taken 
into account. The contribution of brcmsstrahlung to the 
data amounts to 33%. The uncertainty in the correction 
due to internal (external) bremsstrahlung was estimated 
to 2% (1%). Summing these uncertainties quadratically 
leads to a model-dependent uncertainty of 4.7% (3.6%) 
for the analysis using the IPSM (CBF) spectral function. 

The measured nuclear transparency Ta (solid symbols) 
for 12 C is shown in Fig. 0]as a function of the kinetic en- 
ergy of the proton together with previous results obtained 
at SLAC [221 (circles) and Jlab [H, |23| (squares and dia- 
monds). The error bars shown in the figure contain the 
statistical and systematic uncertainty but not the model- 
dependent error. This applies also to the data points of 
the previous works. Since the previous experiments were 
analyzed using the same assumption and ingredients the 
model-dependent error is the same for them, while it is 
somewhat lower in the case of using the CBF spectral 
function. 

The solid lines drawn in Fig. 0] are the result of the 
theory presented in this paper. For comparision also re- 
sults from previous experiments 0, 122I l23l | for iron and 
gold are shown. For all three nuclei and large proton ki- 
netic energy (> 1.5 GeV) the theory describes the data 
well within the error bars. At low energy there is remark- 
able agreement between theory and the experimental re- 
sults obtained using the CBF spectral function. The two 
data points at the lowest T p for 12 C could indicate a de- 
viation from the prediction, but considering the model- 
dependent error bar no firm conclusion can be drawn. 
With the standard analysis the experimental results are 
f=a 5% too low but in agreement with previous analyses 
using the same ingredients. On the other hand the data 
points for gold seem to exceed the theory. For these anal- 
yses a correction factor l/e SRC — 0.78 was used [22l[23l|. 
If one would have used the CBF spectral function the 
results would be lowered by w 7% and thus closer to the 
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theory. 
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